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We offer a unified approach to several phenomena related to the electromagnetic vacuum of a 
complex medium made of point electric dipoles. To this aim, we apply the linear response theory to 
the computation of the polarization field propagator and study the spectrum of vacuum fluctuations. 
The physical distinction among the local density of states which enter the spectra of light propa- 
gation, total dipole emission, coherent emission, total vacuum energy and Schwinger-bulk energy is 
made clear. Analytical expressions for the spectrum of dipole emission and for the vacuum energy are 
derived. Their respective relations with the spectrum of external light and with the Schwinger-bulk 
energy are found. The light spectrum and the Schwinger-bulk energy are determined by the Dyson 
propagator. The emission spectrum and the total vacuum energy are determined by the polarization 
propagator. An exact relationship of proportionality between both propagators is found in terms 
of local field factors. A study of the nature of stimulated emission from a single dipole is carried 
out. Regarding coherent emission, it contains two components. A direct one which is transferred 
radiatively and directly from the emitter into the medium and whose spectrum is that of external 
light. And an indirect one which is radiated by induced dipoles. The induction is mediated by one 
(and only one) local field factor. Regarding the vacuum energy, we find that in addition to the 
Schwinger-bulk energy the vacuum energy of an effective medium contains local field contributions 
proportional to the resonant frequency and to the spectral line- width. 



I. INTRODUCTION 

Dipole emission and the Casimir effects in dielectric 
media are long-standing research topics in quantum elec- 
trodynamics. They constitute paradigmatic phenomena 
of the interaction matter-radiation related to the quan- 
tum nature of the vacuum. Regarding dipole emission, 
it is known since the work of Purcell [lj that the spon- 
taneous emission rate of an atom, T, in a dielectric 
host medium depends on the interaction of the emit- 
ter with the material environment. This is so because 
the back-reaction of the medium on the atom via the 
self-polarization field modifies its self-energy. In the first 
place the host medium modifies the density of channels 
into which the atom can radiate -i.e. the Local Density 
of States (LDOS)- and hence the value of F. Second, 
the integration of the associated self-energy gives rise to 
a shift in the resonant frequency of the atom. Third, the 
dipole-transition-amplitude gets also modified. As a re- 
sult, the medium is said to renormalize the polarizability 
of the emitter. 

In addition and complementarily, the medium polar- 
izes the EM vacuum since material degrees of freedom 
couple to the EM fluctuations. As a result, polarization 
fluctuations arise. In particular, this reflects in a renor- 
malization of the susceptibility of the medium, which in 
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turn determines the vacuum energy. 

We will postulate two distinguishable EM vacua at- 
tending to the existence of two different spectra of fluc- 
tuations. These are, a light vacuum in which the light 
from an external source propagates through the medium; 
and an emission vacuum in which the radiation emitted 
from a dipole within the medium lives. When the dipole 
constituents of the medium are equivalent to each other, 
the spectra of fluctuations in both vacua are related by 
an exact mathematical expression. It will be shown that 
the spectrum of the light vacuum is only part of the spec- 
trum of the emission vacuum. 

In accordance to the above distinction, pairs of ho- 
mologous quantities ascribe to each vacuum. These are, 
the light propagation spectrum and the dipole emission 
spectrum; the bulk-Dyson field and the dipole field; the 
Schwinger-bulk energy and the total vacuum energy. The 
relationships between the members of each pair are de- 
termined by local field factors. The finding of analytical 
expressions for these relations is the main achievement of 
the present work. 

We will deal with a dielectric medium made of isotropic 
point dipoles. The distribution of dipoles is statistically 
homogeneous and isotropic. Ours is a classical Green's 
function based approach. The system is linear, being 
characterized by a transference matrix T. The cou- 
pling matter-radiation is weak so that in any multiple- 
scattering process photons are asymptotically free in be- 
tween any two scattering events. Therefore, the propa- 
gator of the radiation in the medium can be computed in 
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terms of a stochastic Lippmann-Schwinger equation over 
the ensemble of scatterer configurations. Our approach 
is compatible with the existence of internal resonances in 
single scatterers and geometrical resonances in clusters. 
However, excitons proper of high order systems in the 
strong coupling regime [1] are disregarded. 

The medium under consideration is stationary with re- 
spect to (w.r.t.) the emission processes. All the quan- 
tities computed are time-independent ensemble-averaged 
quantities. No other mechanism but radiation is con- 
sidered in the processes of gain and loss of energy by 
the dipoles. In particular, collisional effects and the 
Doppler shift are disregarded. For the latter, the condi- 
tion v/c <C V /uj is required, v being the typical velocity 
of the scatterers and oj the frequency of interest. There 
is also the weak-coupling condition w 3> T. Finally, in 
order to treat the scatterers as point dipoles we demand 
that ui/c <C a -1 and £ > 3a, a being the typical radius 
of a dipole and £ being the minimum distance between 
them. 

Our approach is based on the microscopical diagram- 
matic treatments carried out by Foldy [~3] , Lax Q , Frisch 
Q, Dzyaloshinskii, Gorkov, Abrikosov and Lifshitz [6, 7], 
Bullough and Hynne [1, Q. We base on Milonni's work 
[l"ol - [l"2j to apply our classical Green's function formalism 
to the computation of quantum quantities. 

The article is organized as follows. In Section |TT] 
we derive the fundamental relation between the Dyson 
propagator and the polarization propagator. The self- 
polarization propagator determines the radiative correc- 
tions which renormalize the single-particle polarizability. 
This issue is addressed in Section IIIII together with the 
derivation of an expression for the optical theorem in 
complex media. That expression yields the total dipole 
emission, which is decomposed into coherent and extin- 
guished components in Section [IVJ In Section [V] the 
light vacuum and the emission vacuum are related via 
symmetry transformations. The role of the local field 
factors in coherent emission and field renormalization is 
explained. In Section IVII we calculate the vacuum en- 
ergy of a generic complex medium using the Schwinger 
variational method and particularize to the case of an 
effective medium. In Section IVIII we comment on the 
difference between LDOSes and discuss some misconcep- 
tions present in previous works. In Section [Villi we offer 
an alternative derivation of the fundamental relation be- 
tween Dyson's and the polarization propagators based on 
the Coupled-Dipole-Method. The most relevant results 
are summarized in the Conclusions. 

Regarding notation, we will label three-spatial- 
componcnt vectors with arrows and three-by-three ten- 
sors with overlines. We will denote the Fourier-transform 
of functionals with g-dependent arguments instead of the 
r-dependent arguments of their position-space represen- 
tation. 



II. THE LOCAL DENSITY OF STATES 

In application of the fluctuation-dissipation theorem, 
the LDOS at a point r for photons of frequency oj and 
any polarization is given by the imaginary part of the 
trace (Tr) of the w-mode of the propa gato r for virtual 
photons created and annihilated at f [71 Il3l - fl6| , 



Af(r» = -^{Tr{G(f,f;u/)}}. 



(1) 



When dipoles are present in the medium, in addition to 
the in-free-space EM fluctuations, polarization fluctua- 
tions arise as a result of the interaction matter-radiation. 
Consequently, the EM vacuum gets polarized and the di- 
electric matter acquires additional self-energy. The field 
whose fluctuations in vacuum are named EM vacuum 
fluctuations is the s elf-p olarization field, also called radi- 
ation reaction field [lCf. 

Amongst the polarization fluctuations, polaritons are 
excitations which satisfy certain dispersion relations and 



propagate coherently 17J. That is, they are normal 



modes and they determine the spectrum of coherent 
dipole emission. We will show in Section IIVI that coher- 
ent emission consists of two components. Those are, a 
direct component which is transferred radiatively and di- 
rectly from an emitter into the medium. And an indirect 
component which corresponds to the coherent emission 
radiated by the induced dipoles surrounding the emitter. 

We will denote the total LDOS accessible to the pho- 
tons of frequency oj emitted from a dipole at r by 
Af emzs ' (r,uj). We will show that the aforementioned di- 
rect coherent emission has the same spectrum as the ra- 
diation which propagates through the medium from an 
external (uncorrelated) source. Thus, we will denote 
by N h3ht (r, oj) that part of M emls - (r, oj) accessible to 
monochromatic external light of frequency oj. As long 
as all the dipoles have the same polarizability, any point 
f is statistically equivalent to any other so that we can 
drop the r-dependence on both LDOSes in this case. 



The propagator of the coherent-Dyson field and 

flight 



A generic random medium is made of N + I point 
dipoles or scatterers in a volume V such that, in the 
thermodynamical limit A, V — >■ oo, the average numeri- 
cal density p = is finite. Each dipole is labeled by 
a subscript i which takes integer values in [0, A]. Each 
scatterer configuration is determined by a set of A + 1 
scatterer position vectors {Ri} and a set of A + 1 renor- 
malized isotropic polarizabilities {a 1 } such that a 1 = a l I 
-see below for the definition of renormalized polarizabil- 
ity. In a statistically homogeneous and isotropic medium 
all the dipoles have a common polarizability, a 1 = a 
Vz = 0, ..,A, and the position vectors are statistically 
equivalent as performing ensemble averages. We will 
stick to a statistically homogenous and isotropic medium 
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unless otherwise specified. 

A/" 9 is given by the Green function of the ra- 
diation which propagates from an external source (i.e. 
ideally located at infinity) through the medium to any 
point. It is sketched as a series of multiple-scattering 
diagrams in FigHJ Because both the external source 



and the end point are spatially uncorrelated w.r.t. any 
of the scattering events in all the diagrams, the corre- 
sponding Green's function or propagator is the radiative 
component of Dyson's, G(r,r';uj) [g, @]- The stochastic 
Lippmann-Schwinger equation for G(f,r';uj) reads @, 



G ij (r,f;uj) = G i ° ) (f-f;uj) + 



V dV" G$ (r - r" ; uj) (t% m (/ , ?' ) ) G mj (/ 



(2) 



where G (0) (r - r 1 ; uj) is the in-free-space propagator, 
k = uj/c and summation over repeated indices is im- 
plicit. The functional t u (r , r ) denotes the t- matrix 
of each specific configuration. It relates to the auxiliary 
matrix i u (Ri,Rj) through 



t u {r ,r ) 



N 

E 

i,j=0 



t w (Ri,Rj) <5 (3) (r" - Ri)S^(r" - Rj). 

t"(Ri,Rj) propagates photons along all possible 
multiple-scattering trajectories connecting the dipoles i 
and j at points Ri and Rj. It has a cluster expansion, 
l"{Ri,Rj) = -k 2 &5 l3 + E„= 2 ^ n) (^>4-Wi 5 where the 



Ah 



term in the sum contains n factors — k 2 a and n — 1 



tensors G^°>{R m — Ri;uj). As an example, 

N 

tW(iMj) = (~k 2 &) 4 Y, G^(R\ - Rf,ui) 

■ G^ a \R\ — R m ;uj) ■ G<® (R m - R 3 ; uj) . (4) 

The big brackets in Eq.((3|) stand for the average per- 
formed over the ensemble of scatterer-configurations. 
The subscript 1PI signals the restriction to one-particle- 
irreducible (1PI) correlations. Thus, the stochastic kernel 
of Eq. is the electrical susceptibility function made of 
the sum of 1PI multiple-scattering processes 0, H| , 



-k~ 



t u {r ,r ) 



ipi 



(5) 



The average process may be performed over both the 
classical and the quantum degrees of freedom of the con- 
figurations. In the simplest case, only spatial correlations 
among the scatterers are considered Q. 

The in-free-space propagator, G^ (r — r 1 ; uj) : is the ex- 
pectation value of the electric field generated at the point 
r by a classical dipole at r 1 oscillating with frequency uj. 
Equivalently it is the w— mode of the Green function of 
Maxwell's equation in free space, 



— r"; uj) has also a quantum-mechanical interpreta- 
tion. It is the scattering amplitude, computed at second- 
order of perturbation theory, of the interaction between 
two point dipoles of unit dipole moment placed at r and 
r* respectively -see eg. [ljj. G^\r;uj) consists of an 
electrostatic (Coulombian) dipole field propagator, 



-1 

47r r 



and a radiation field propagator. 

^ikr _ ^ 



a i kr 



—4irr 



72 V ® V J A • 



(7) 



(8) 



In the reciprocal space and for an isotropic medium, any 
tensor can be decomposed into longitudinal and trans- 
verse components with respect to the propagation direc- 
tion, q. In free space, 

m (q; uj) = G { 1 ] (q; uj) (I - q ® q) + g[, 0) (q; oj)q®q, 

where q is the unitary vector parallel to q and 



Gf(q;u>) 



k 2 — q 2 



Gf\q;uj) = 



k 2 ' 



(9) 



While the radiative component is fully transverse, the 
electrostatic one is fully longitudinal. 

For a given medium, G(r, f 1 ; uj) is the expectation value 
of the coherent-Dyson field or bulk field, Ep(r), gener- 
ated at r by a source at f", being both points uncorrelated 
w.r.t. any of the scatterers in the medium. In particular, 
when r* is located at infinity, it is the transverse com- 
ponent w.r.t. the propagation, G±[f, f"; uj), that propa- 
gates coherent external light [8| . G is obtained averaging 
over the ensemble of scatterer configurations all possible 
multiple-scattering propagation trajectories which start 
at r 1 and end at r. G(r,r t ;uj) is said to be the analog of 
G^°'(f — r'^uj) in a complex medium in the sense that it 
satisfies an ensemble-averaged Maxwell equation analo- 
gous to that in Eq.©, 



V x Vx 



G(°'(r-r';u) = (f (3 '(r , -r')I. (6) / 



-e"(r)l- V x Vx 



gir^-uj)^ = (10) 
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In the above expression e u (r) = 1 + J2iLo o^ 3 ' (? ~ 
Ri) and g(r,f';uj) are respectively the dielectric ten- 
sor and the Green function of each specific scatterer- 
configuration. The integral version of Eq. (fTU|) is that of 

Eq.© for G{f,f\u) = (g(r,f*;u))\. 

In a statistically homogeneous medium G(r, r"; uj) and 



X u '(r,r') are functions of 



for any pair of points. 



In the Fourier space, isotropy allows to split the Dyson 
equation for G(q) in two uncoupled and mutually orthog- 
onal algebraic equations, 



G Ul (q) = G^ l ( q ) 



k 2 G 



?},(?) Xx,||(?)G Xl ||(«). (11) 



In the above equation and in the following we will omit 
for brevity the explicit dependence of the functionals on 
to. The same as for the i-matrix, X-LJI (?) adjust to cluster 
expansions of the form, 



X±,\\ 



(12) 



n=l 



The functions -X"j™|i (?) incorporate the spatial dispersion 
due to the 1PI spatial correlations within clusters of n 
scatterers. In addition, they account for the intermedi- 
ate multiple-recurrent-scattering processes among them. 
In field theory terminology, x is proportional to the self- 
energy tensor, £ji|(<z) = — k 2 \±.\\ (<?)• Alternatively, 
Eq. (1111) can be written in terms of the ensemble-averaged 
T-matrix, 

G±,\l(q) = Gf Al (q) + G { °\( q ) T ±i]l (q) G { °\(q), (13) 

where, in application of isotropy and homogeneity, Tjii 
are defined through 

t%?)) = (2n) 3 6^(q-q')f(q), 

f(q) = T x (q)(l-q®q) + T {[ (q)q®q. 

The solutions to Eq. (fTTj) are the Dyson propagator com- 
ponents, 



G±(q) = 
G\\(q) = 



1 



k 2 [i + xAq)]-q 2 ' 
l 

k 2 [l + X\\(q)Y 



(14) 



In terms of free propagators and self-energy tensors, 
G can be depicted perturbatively as in FigJTJ Normal 
modes are given by the poles of G±n(q). They satisfy 
the dispersion relations, 



k 2 e ± (q)-q 2 \ 



e\\(q)\ q =k^- 



0. 

= 0. 



(15) 
(16) 



with ej_ ,|i (9) = 1 + X±.\\(q) being the dielectric functions 
for transverse and longitudinal modes respectively. Lon- 
gitudinal modes need of coupling to matter to propa- 
gate while external light excites only transverse modes 
8]. Thus, we have 



Auj 



M liah \uj) = j 



d 3 g 

(M 



rG 



.(?;")}• (17) 



= G I'- 



FIG. 1: Diagrammatic representation of the Dyson propaga- 
tor, G. 



B. The propagator of the polarization field and 



The Green function which determines jV emls ■ (uj) 
through Eq.([T]) is that of the polarization field, Q. For 
a given medium, G(r,f';uj) is the expectation value of 
the electric field generated at the point r by a classical 
dipole oscillating at frequency uj at a point f 1 within the 
medium. It is in this sense that Q(r,f , ;u>) is the ana- 
log of G^°'(r — f'fUj) in a complex medium. Because, 
by assumption, our medium is made of indistinguishable 
point scatterers we choose the dipole emitter to be the th 
scatterer out of the N + 1. Thus, the ensemble-averaged 
Maxwell equation reads here, 



T e u (r)l-VxVx g^r'jw) 



_ = S^(f-R )l. 
=Ro 

(18) 

The restriction | f »_^ signals the fact that the ensemble- 
averaged is performed keeping r* fixed at the position 
vector of the th scatterer. This simple but crucial 
fact makes Eq.([Tg|) different to Ea.([TD]). The resul- 
tant integral Lippmann-Schwinger stochastic equation 

for Q(?y\u) = (g{f,f;u}) i 
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d 3 r 



'd 3 r 



1PI 



=Ho 



(19) 



The expression 



-v -> 

r , r 



IP/ 



f=ilo 



(20) 



is a symbolic manner of writing the stochastic kernel. 
The explicit dependence of i w on f* and the restriction 
denote that while performing the average over all 

the 1PI processes, f* in Q(r , f*;u>) is kept fixed at the 
center of the th dipole. This is a consequence of the fact 
that the stochastic operator on the l.h.s. of Eq. ([T8l) is 
correlated to the source term on the r.h.s of that equa- 
tion. In Eq. (fT9|) the dipole at r* acts as a source and 
hence it is a dipole emitter. In general, the presence of 
an emitter would break statistical homogeneity in a ran- 
dom medium and Q would depend explicitly on the emit- 
ter location, R . However in a statistically homogeneous 
medium, since all the position vectors are equivalent as 
performing ensemble averages, we can drop the explicit 
dependence on R in favor of a generic argument r* as in 
Eq. (IT91 . Note that if we replace r* by r in the stochastic 
kernel and fix r at a dipole, that dipole would be a par- 
ticle being polarized instead. 

Our interest is in the self-polarization propagator that 
enters Af emls - (uj). It is computed out of Eq. ([TO|) by set- 
ting It* — f] < a. This way, the dipole emitter is at the 
same time the dipole being polarized. In turn, the dis- 
tinction between r and r* is just formal as it cannot be 
resolved for point dipoles. We will set r — r* without 
loss of generality. Q(f,r;u)) carries the radiative correc- 
tions which enter the renormalized polarizability of the 
emitter, a. Because both the fluctuations of the classical 
self-polarization field and those of the quantum vacuum 
field are related by the fluctuation-dissipation theorem 
-see eg.[l^,[ll|, both classical and quantum corrections 
can be treated within the same formalism. Clarifications 
on this point will be provided where necessary. 

The explicit dependence of Eq. (|2"0)) on the emitter po- 
sition in Eq. p^)) makes that any scattering process of 
the virtual photons in their way from and towards the 
emitter be correlated w.r.t. the emitter position. When 
the emitter is in all equivalent to the rest of host scat- 
terers, such a correlation is the same as that among the 
host scatterers themselves. In the simplest case the only 
correlation is that of an exclusion volume. It is in this 
sense that the homogeneity of the medium is only virtu- 
ally broken and the exclusion volume around the emitter 
is referred to as virtual cavity. Should the emitter be 
distinguishable w.r.t. the other scatterers, that correla- 
tion would differ. Homogeneity would be actually bro- 
ken and the cavity would not be virtual but a real cavity. 
The breaking of homogeneity, either virtually or actually, 



makes the computation of Q(r,r';uj) different to that of 
the ordinary Dyson's propagator and the stochastic ker- 
nel of Eq. p^l different to the susceptibility tensor. We 
will refer generically to the scenario in which the emitter 
is in all equivalent to the host dipoles as stricto sensu 
(s.s.) virtual cavity (VC) scenario. It is for this scenario 
that there exists a close relation between G(r, r; uS) and 
Q{r, r;oj) which we proceed to investigate. 

In the following, we will work in Fourier space and 
write G(r, r;ui) as 



G(r,r;u) 



d q n , 



(2tt) 3 ' 

[2^ (0) + v^ (0) ] + + <pT c ] I ( 21 ) 



where yj| are the in-free-space values which contain 
real divergences to be regularized. ^""m are the di- 
vergenceless scattering pieces. Physically, the ^-factors 
account for the dipole self-energy associated to self- 
polarization photons. The numerical prefactor 2 in front 
of ip±_ stands for the two transverse polarization modes. 
The w-dependence will be omitted hereafter unless nec- 
essary. We first observe that Q{r,r) is made of 1PI dia- 
grams in which the end points coincide. Those diagrams 
amount to the so-called recurrent scattering. Let us draw 
them as linear processes as in Fig(5] The second observa- 
tion is that, by reciprocity, every intermediate scattering 
event in Q (r , r) which is correlated to the fixed dipole at 
r will be correlated either at the near end or at the rare 
end of each diagram, indistinguishably. Taking advan- 
tage of this feature in every multiple-scattering diagram 
like that in Figj2]J6), we can attribute all the irreducible 
correlations of the intermediate scattering events to the 
emitter on the left. By proceeding so, we end up with 
an effective separation of all those pieces irreducibly cor- 
related to the emitter on the left completely untangled 
from those non-lPI pieces on the right. The sum of the 
1PI pieces on the left amounts to x/p^: where the factor 
1 1 pa, stands for the removal of the first random scatter 
which enters the diagrams of x m favor of the fixed emit- 
ter location. The sum of the non-lPI pieces on the right 
amounts to the bulk propagator G. As a result, we end 
up with the formulae -FigsHJc, d), 



Ql c {i) 



G ± (q) 



x±(q)/(pa) 
k 2 [i + xAq)}- 



pa 

xn(g) r (n) = — *ii (g) 
pa 11 W pa fe2[l + X ||(?)] 



j 7 (22) 
(23) 
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■■h(r) ! \=8 { \) %= pn 



(6) 



• (c) 



(d) 



In terms of the <p- factors defined in Eq. (|2"Tj) . the LDOS 
of emission reads, 

M emls -{u:) = -|^3|2^ (0) + 2i P u A sc + ipf c }. (30) 

Note that, differently to the case of AT b9ht (uj), even for 
the case that longitudinal normal modes do not exist, 
both transverse and longitudinal modes contribute gen- 
erally to J\f emts - (a;) as it does not solely include coherent 
modes. This being so regardless of the absorptive prop- 
erties of the medium. 



FIG. 2: (a) Feynman rules. Only two-point irreducible cor- 
relation functions, h(r), have been used for the sake of sim- 
plicity. (6) Diagrammatic representation of the equivalence 
between multiple-scattering processes amounting to G(f,f). 
Reciprocity applies, p represents the virtual presence of the 
emitter. (c),(d) Diagrammatic representations of the integra- 
tion in q of Eq, (122p and Eq. (|23[) respectively. 



The above expressions contain both G^ and G^ which 
carry divergences. The scattering terms, 

vc f d3 9 f Xx(q)/{p&) 



(2ir)nk*[l+ X ±(q)]-q 2 



-Gf(q) ,(24) 



d 3 g 



1 



x\\(q) 



pa k 2 [l +X\\(q)] 



-G§\q) , (25) 



are however fully convergent. 

Alternatively, we can write Q (q) in other forms mak- 
ing use of Dyson's equation. Using Ea. (ITT1) we get, 



k 2 pa 



G i 



G 



(o) 
-J-,11 



(26) 



In function of the T-matrix of Eq. p^)) it reads, 



— I 



(27) 



The reader can find in Section I Villi an alternative proof 
of Eq. (|27|) using the formalism of the Coupled-Dipole- 
Method. 

Finally, we write the Lippmann-Schwinger equation of 
Eq. (II9l) for the polarization field propagator in Fourier 
space, 

Ql%q) = G% (q) + Gf\ (q) E^(q) (28) 
where 



-pa 



l w r>(°) 
X_l,|| G v 



J-,11 



X±,\\ 
pa 



+ fc a Xx,||Gi°J,J(5) (29) 



is the stochastic kernel of Eq. (f2"0")) expressed as a function 
of Xi G^ and pa. 



ILL RENORMALIZATION OF THE 
SINGLE-PARTICLE POL ARIZ ABILITY 

The propagator Q{f, f)oS) carries the radiative correc- 
tions which renormalize the so-called bare-electrostatic 
polarizability ag of a given dipole. «o can take different 
forms depending on the nature of the dipole. We refer 
to [III for a comprehensive summary. In Appendix [X] we 
adopt a semiclassical model in which dipoles are spherical 
nanoparticles of permittivity e e (u>) and radius a <C c/cj 
with ao(uj) = 47ra 3 ^ e [^] + 2 ■ in Appendix iBl we adopt a 
two-level atom model in which the atom is in a spheri- 
cally symmetric quantum state. 

Let us consider first the classical model and compute 
the off-resonant renormalized polarizability of a classical 
dipole. Let us suppose that the dipole is stimulated by 
a stationary external field which oscillates in time with 
frequency u>, Eq. The processes of absorbtion and emis- 
sion of radiation by the emitter become stationary after 
a relaxation time of the order of In function of the 
95-factors, the average power emitted and absorbed by 
the stimulated dipole reads -see Appendix lAl 



W, 



2 l I + i 



fc 2 a [2^ 0) +2^ c + ^ 



-}\E \ 



-w 3 e 



|ao| 2 3{2^ ( L ° ) + 2^ + ^} 



6c 2 



1 1 + \k 2 a [2^l } + 2<p s f + ipf ]| 2 



(31) 



%{a } 



fc2 \l + \k 2 a [2<p 



(o) 



sc] 1 2 



l^o| 2 (32) 



The explicit dependence of the ip- factors on k = lo/c has 
been omitted. The term in Eg. (|3"2"]l corresponds to the 
power absorbed within the emitter while that of Eq. (|3"Tj) 
corresponds to the power transferred to the medium. 

We can write W^ ot in terms of a renormalized polar- 
izability a as Wl ot = ^%{aE Q ■ i%} = af-\E \ 2 $s{a}, 
with 



a(k) = 



1 + |fc 2 a ^i 0) + \k 2 a [2y*£ + 



(33) 



ip- 



Because the computation has been assumed for frequen- 
cies far from resonances, it is possible to regularize the 
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real divergence of ipy by setting it to zero so that 

2ifY — m the above equation. Finally, writing 

Eqs. (|31l32j) in the form, 



w ^ ot = -^i^i 2 [^ { ? + + ^ c } 



3 3{a } 



a 



by means of the fluctuation-dissipation theorem. While 
it can be interpreted that only the latter contributes 
to Eq. ([34]) . both of them, and in the same fashion, 
contribute to Eq. (j56"|) . Following Milonni's arguments 
, JIol [IH, this explains why the Einstein coefficient for 
" spontaneous emission is twice that for stimulated emis- 



we obtain an expression for the optical theorem in com- 
plex media. 

Next, let us consider the polarizability of a quantum 
two-level atom. In good approximation, we can adjust a 
to a renormalized Lorentzian polarizability [20j , 



IV. ON THE NATURE OF STIMULATED 
EMISSION 



f 



a = -&„k res [k res -k -iTk /(ck res )] 



(35) 



This way, a gets parametrized in terms of the renormal- 
ized values of <So, k res and T. Radiative corrections enter 
in the same fashion as for the classical model above 
-see also Appendix [Bj Therefore, the parametrization in 
Eq.p5p allows us to give a physical meaning to the tp- 
factors and to their in-free-space divergent terms. In free 
space ao relates to the dipolc transition matrix element 
p and the bare resonant frequency, u>q = cfco, through 

2 1 1 2 

Also, the above parametrization yields 



an automatic regularization of the divergent real part of 
2ip ^ , $i{2(p^} = • This regularization procedure 

determines the in-free-space Lamb-shift when coupling 
to bare radiative modes is considered [21] . By consis- 
tency, To = caofco/67r = 3 ^° ah \p\ 2 is the in- free-space 
decay rate. Comparison of Eq. (l33l) with Eq. (|35|) yields 
the following expressions for the renormalized parameters 
in terms of the <p- factors, 



r = -^ fc 3 3{2^ 0) 



2tt 

'T2' 



(0) 



2cpf 
-2^1 



(36) 



where k res is a real non-negative root of the equation 



(fc/fco) 



1 



and a = a (k /k res y 



(37) 



(38) 



Except for Eq.(|38|). the formulae for k res and T have been 
already derived using a fully quantum-mechanical (QM) 
formalism [2~iT - [2~4l | . In [23|, [24| the same self-consistency 
requirement between Eq. dMl) and Eq. (|57j) was obtained. 
Note however that only the resonant term of the level 
shift is accounted for in Eq. (|3"7| [231 ] . That is, we have 
implicitly assumed that the level shift due to the van-der- 
Waals interactions is the same for the two atomic levels 
[13]. See also Section ED 

The underlying reason why the radiative corrections 
which yield the spectrum of stimulated emission in 
Eq. dMl) enter in the same fashion the spontaneous de- 
cay rate in Eq. dMl) is that the spectrum of the self- 
polarization field relates to that of the vacuum field 



Let us consider the emission term of Eq. (pH) for a sce- 
nario in which only one of the classical dipoles of a sta- 
tistically homogeneous dielectric is stimulated by an ex- 
ternal field. We will refer to the emitter as dipole source 
as well. The (/3-factors are those for the virtual cavity 
scenario and we can write, 



Wl ot = W 



d 3 
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(2nY 



Tr- 



( 9{ ^M.G(g; W )}}, (39) 
(pa J 



where W Q 



6c eo 



\po\ , being po = e^aE^r) the dipole 



moment induced by the monochromatic external field 
Eo(r) of frequency us on the emitter at r. 

The first obvious decomposition is that between trans- 
verse, Wj°*, and longitudinal emission, W{; ot , 



W, = 



W 

Wo 



/|2 3{ ^ Gi ( 9 ) } , (40) 
[Ziry pa 

d 3 ? X|(qO 
——^^{^-^G^q)}. (41 
(27r) J pa 



Hereafter we drop the script w bearing in mind that all 
the quantities depend on the frequency of the external 
field. 

Coherent emission is that carried by the coherent- 
Dyson field itself and by any other averaged-field com- 
ponent in-phase with it. That is, the power carried by 
the total ensemble-averaged field radiated [3]. This is 
the radiation whose modes satisfy the same dispersion 
relations as the normal modes in the bulk, Eqs. (|15ll6[) 
They are the poles of G±(q) and G^q) in the 
integrands of Eqs. (|4QI41[) . The rest of the radiation goes 
into dispersion and absorbtion. It is termed generically 
extinguished emission [26| . Therefore, we define, 



W Coh. 



wr 



Coh. 



IF, 



W\\ xt - 



Wo 
Wo 
Wo 
Wo 



<] q 25R{^M}3{G ± (g)}, (42) 



(27T) 



pa 



|« 29{ ^M }R{Gl(?)} , (43) 
(27r) J pa 

sft{^}Z{G „(«)}, (44) 
(27r) J pa 

|L a{ MM }Sf { G|| ( 9)} . (45) 
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The definitions of formulae Egs. (|42l44[) as trans- 
verse/longitudinal coherent emission are inspired by the 
spectrum of the on-shell modes of field theories p], [27| ■ 
Hence, the sourceless counterpart of W^ oh - = W±° h ' + 
B|l ■ is the spectrum of polaritons. In the following, we 

examine W^ oh " in the framework of Classical Optics to 
show that it corresponds indeed to the ordinary classical 
interpretation of coherent radiation. Let us write W% 



Coh. 



6c 2 eq 



(46) 



r'Tr{K{x(F- r)/pa} ■ %{G(r*,?)}}, 



where both \ an d G are written in the position-space rep- 
resentation for convenience. Considering the fields classi- 
cally, the fluctuation-dissipation relation reads @, EH, Ell , 



9f{G(r 



;)} = 



hk 2 



(EZWQEgtf), (47) 



where Ep(f) is the w-mode of the 
and the script D stands both for 
emission for the reasons explained 
and writing ^R{x{r — r 1 )/ pa} = 
5(V(r~7)t} + S^(f- f)I in Eq 
plicitly the field radiated directly 
from that which is emitted by the 
dipoles, 



coherent-Dyson field 1 
Dyson and for direct 
below. Using Eq.([47]l 
[3?{x(f - f)/pa} - 
PB^l . we separate ex- 
by the source dipole 
induced surrounding 



Coh. 



TTUJ 



^-\p Q \ 2 (E" D (r) -E%*(r)) 
^|p | 2 Tr{|dV([5R{x(r'-r')/pa} 
S^(r- ■ ® E^(f))}, (48) 



where the brackets here denote vacuum expectation val- 
ues. The first term on the r.h.s. of Eq. (|48| contains 
the coherent power carried by the field directly radiated 
by the dipole source into the medium 28], W^ oh " — 
^\po\ 2 (\Ep(r)\ 2 ) . transverse component is the same 
expression as that for the coherent intensity carried by 
a propagating light beam whose source sits out of the 
medium In the rest of the r.h.s. we can identify part 
of the field emitted at f 1 by the induced dipoles sitting 
around the dipole source which propagates towards the 
source itself located at r, 



E?{r)= d 3 r'[x(r-f f )/pa-5^(f-f f )I]-E%(r'). (49) 



1 Throughout this paper the electric field vectors with script u) 
are frequency modes of the field, not to be confused with the 
monochromatic external electric field of frequency u) used in Sec- 
tions [rrrjviii] 



The subscript / stands both for induced and for indirect. 
Therefore, we can write 



w: 



Coh. 



= g^ol 2 (\E%(r)\ 2 ) +mE?(r) ■%?(?))} 



and W; 



Coh. 



(50) 



W Coh. +W Coh.^ with 



W, 



Coh. 



= —\p \ 2 ^{(Ef(r)-E^(r))}- (51) 



As expected, Eq . ([ST]) is the coherent power of that com- 
ponent of the averaged-field induced on the surrounding 
dipoles which is in-phase with Dyson's -see [29| for an 
analogous computation in a simpler scenario. 
Alternatively we can write, 



w: 



Coh. 



D 



Coh. 



6c 2 eo 

-UJ 3 

6c 2 eo 



\p \ 2 Tr{<Z{G(7,7)}}, 



(52) 



Wo 



-'Tr|3?{x(r - f)/pa 



- 6W(r-T*)I]}-S{G(r',r)}\. (53) 

In the diagrammatic representation of W£ , G is di- 
rectly attached to the emitter. 3{G^} yields the trans- 
verse normal modes or actual photons which mediate 
the radiative energy transfer from the emitter into the 
medium [18j. The diagrammatic representation of Wf oh " 
can be split in two pieces. The induction piece, x(r — 
f 1 )/ pa — J*- 3 ' (r — f")I, is attached to the emitter. Its non- 
trivial real part, ffi{x(r — r')/pa}, does not possess any 
normal transverse mode. In particular, no actual pho- 
ton emerges from the emitter, being the induction due to 
non-radiative energy transfer mediated by virtual pho- 
tons. The second piece amounts to G and is attached 
to the host dipoles. 3{Gj_} yields the transverse normal 
modes or actual photons emitted by the induced dipoles 
which are in phase with those of W£ oh ' . We conclude 
that only the spectrum of the power transferred radia- 
tively and directly into the medium is equivalent to the 
spectrum of external light. 

In the extinguished emission the factors Sj X± -V q ' > } 
of Eqs. (l43l45[) contain both absorbtion and dispersion, 
which are in principle indistinguishable observationally. 
The former is attributed to the intrinsic imaginary part 
of the bare polarizabilities, 3{ao}- The latter is captured 
by the bare transverse modes affected by correlations 0] . 
That is incoherent radiation (in the classical sense) which 
cannot be written in a form analogous to that of Eq. (I5T1) . 

Observationally, the power collected in the far field at 
a distance r' £ from the emitter -£ being the typical 
correlation length- using an integrating sphere would be 
the sum of the coherent and the dispersed emissions. On 
the contrary, the coherent component would be propor- 
tional to the integral of the square of the averaged-field 
measured by an observer on a sphere of radius r' . In an 
effective medium it relates to WY° h ' through the Beer- 



9 



Lambert attenuation factor 2 , 

W Coh -{r') = W^-expi^nujr'/c} 



cneor 



.'2 



E(r',Q)) dfi, (54) 



where n, k are the refractive index and the extinction 
coefficient of the medium respectively and the big angular 
brackets denote the ensemble-average. 



V. RELATION BETWEEN LIGHT 
PROPAGATION, DIPOLE EMISSION AND 
THEIR VACUA 

A. The vacuum of light vs. the vacuum of emission 

It is plain from the above equations that only in free 
space G(r,r) and G(r,f) are equal to &°\r,r) and so 

are N Ught (u) and Af emls -(oj) to Af°(u) = On the 

other hand, the fluctuation-dissipation theorem relates 
the quadratic vacuum fluctuations of the electric field 
operator at a point r with the imaginary part of its prop- 
agator from r to r @, [H, HH , 



hw 2 
e c 2 ir 



3{G-(r, r; 



)}. (55) 



Therefore, we infer that the photons emitted by an iso- 
lated excited dipole propagate in the same vacuum as 
the radiation which propagates from an external source 
in absolute absence of dipoles. This is obvious since the 
source is an isolated dipole itself. However, because in 
a complex medium G(f, r) ^ Q(r,r), we infer that the 
EM vacuum in which external light propagates, \£l) h9ht , 
is different to that of the photons emitted by the point 
dipoles, \ fl) emls -, 



light {VL^^r) <g) El(f)\fl) li9ht = 

emis - (n\E u {f) ® El(r)\n) emis - 



-hu 2 
e c 2 7r 



3{G ± (r,r»}(56) 



-km 2 
e c 2 7r 



Let us consider the s.s. virtual cavity formulae of Sub- 
section QlB] and let us represent polaritons and EM fluc- 
tuations in general by closed loops of virtual photons cre- 
ated and annihilated at the location of a probe dipole, r. 
We can compute the n-scattering loops which amount 
to Q vc {r 1 r) out of (n + 1 )-scattering loops with unde- 
fined reference frame -see Figj3ja). To do so, r must be 
chosen amongst the n + 1 position vectors of the scatter- 
ers and we must treat the scatterer chosen as a virtual- 
probe emitter to be removed. If r is uncorrelated w.r.t. 



2 Note that in Eqs. (I47I51II the electric fields are operators acting on 
the vacuum state. In Eq J54H the electric field is an expectation 
value. The ensemble-average is performed over that expectation 
value. 



any other point along the original loop, the resultant po- 
lariton amounts to G(r, r) and hence, to G vc (r,r) -see 
Figj3j6). That is the same situation met by the photons 
coming from an external light source at infinity. On the 
contrary, if f is correlated to any of the remaining dipoles 
in the original loop, the resultant polariton amounts to 
G vc (r, r) but not to G(f,f) -Fig[3jc). Translation in- 
variance holds statistically in the sense that it is the ori- 
gin of the probe-dipole reference frame, r, that can be 
any point in the medium. 

The finding that the polaritons in G(r, r) are all in- 
cluded in Q vc (r, r) was already obtained through the 
study of coherent emission. Another way to look at this 
relation is by considering the symmetries of each vac- 
uum. If the origin r of one of the n-scattering loops of 
Qvc ^j? ^ j g translated to any of the remaining n scat- 
terer position vectors of the original {n + I)-scattering 
loop, the resultant loop is also a polariton that amounts 
to G (r, r). On the contrary, if the same transforma- 
tion is carried out on a given n-scattering loop of G(r, f), 
there is a chance that the resultant polariton amounts 
to G v (f,r) but not to G(f,r). Only those transforma- 
tions which transport f to another uncorrelated scatterer 
generate a loop in G(r,r). The group of transformations 
which leave G{r, r) invariant is the symmetry group of 
flight and likewise for G vc {r,r) and \Q) emls - . There- 
fore, we conclude that the group of symmetry of the light 
vacuum is included in that of the emission vacuum. 




G 



FIG. 3: (a) Five-scattering closed-loop with undefined refer- 
ence frame, (b) Fixing the reference frame at the position of 
the uncorrelated (removed) scatterer labeled as I yields a po- 
lariton in G(r, r). It can be excited either by external light or 
by dipole emission, (c) Fixing the reference frame at the po- 
sition of the correlated (removed) scatterer labeled as 3 yields 
a polariton in Q vc (r,r) which is not in G(r, r). It can be 
excited by dipole emission only. 
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B. Local Field Factors as renormalization functions 

The relation of proportionality between Dyson's and 
the polarization propagator in Fourier space, 



Gx,\\(q;u>) = Cx,\\(q',u) Gx,\\(q;w), 



(58) 



serves us to define C\ u « (q; u>) 



(9) 



as the local field 



■'± t \\\i7—j p a 

factors (LFFs) of the s.s. virtual cavity scenario. A rela- 
tion of the same kind in the specific scenario of emission 
from an Onsager cavity drilled in a continuous dielec- 
tric was found by Tomas [3(| and Dung et al. [3l[ for 
long-wavelength modes. In the present case the propor- 
tionality relationship holds for the whole g-spectrum and 
the result applies to a generic random medium. We now 
recognize that the induction process which enters Wf 
in Section [IV] is mediated by LFFs. 

We discuss here in field-theory terms the proportion- 
ality relationship between the spectrum of light and that 
of coherent emission, which derives from the expression 

Q^lffou,)} = sft{^^}3{G ± ,||(<?;a,)}. From our 
derivation of the direct and indirect components of the 

' Xx,n(9) ' 



coherent emission we know that 3? 



take ac- 



count of the actual polarization due to the closest scat- 
terers around the emitter. An analogous phenomenon 
takes place with the virtual polarization of the vacuum 
in Quantum Electrodynamics (QED) -see eg.[32j and 
Ch.7,10 of [27]. Further integration of Q{g^f-(q; w)} in 

— «3 f X X,|| v^xTll') 



corre- 



momenta gives the residues = Sft| X± ' n pa ± '^ - j 
sponding to the on-shell modes propagated through the 
bulk by Gj_ \\ (q; ui). Thus, we can write the propagator 
of coherent modes in a more familiar way, 3 



Z^G ± , 



few). 



(59) 



In field theory terminology, the functions are renor- 
malization functions. If we stick to the QED inter- 
pretation, the electromagnetic field can be renormal- 
ized in order to get rid off the {Z p ^) n factors which 
arise otherwise in the vacuum expectation values of the 
products of 2n EM averaged-field operators. In doing 
that, the electric field so renormalized is the Dyson field, 

^lx,|| = TO]- 1/2 ^ oh .U,||. 

It is perhaps worth repeating to avoid confusion that 
the coherent nature of the dipole field here is classical. It 
is only the formalism borrowed from QED that is quan- 
tum. 



The possibility that —4 11 present singularities is disregarded 
here. We do not yet understand their possible physical mean- 
ing. It is our guess that it could signal the transition to strong 
coupling. 



C. Coherent emission in a Maxwell-Garnett 
dielectric 

We illustrate the decomposition of the emission spec- 
trum carried out in Section IIVI with the computation of 
the spectrum of coherent emission in a Maxwell-Garnett 
(MG) dielectric. By definition, an MG dielectric is a 
random medium made of well-separated dipoles whose 
spatial distribution is sufficiently well-described by a two- 
point correlation function of the form 0, [26| , 



.MG 



(r)=-/(r-0 



— 1 for r < £ 
for r > £. 



f(r — £) is this way a spherical-exclusion- volume cor- 
relation function of radius £. Usual forms are those of 
a Lennard- Jones potential and a hard-sphere potential. 
For long wavelength modes, q(; -C 1, and low frequencies, 
!:( « 1, the detailed shape of f(r — £) is irrelevant and 
the quasicrystalline approximation becomes exact in the 
computation of the effective susceptibility (33|. Thus, it 
takes the familiar form, 



pa 



1 — ipa 



(60) 



We omit for brevity the explicit dependence on uj. Beside 
the long wavelength and the low frequency restrictions, 
the MG formula disregards an inherent self-correlation 
term in h(r) which gives rise to recurrent scattering. For 
off-resonant frequencies Eq. (1601) is a good approximation 
up to 0[(pa) 3 ]. Bearing this in mind, we proceed to com- 
pute the LDOS of coherent emission, J\f Coh - , accessible to 
a point emitter in an MG dielectric. M Coh - is the LDOS 
of the observable emission which adjusts to Eq. (l54|) . The 
emitter can be either one of the dipole constituents of the 
medium or an interstitial dipole 20]. Emission must be 
stimulated with frequencies far from internal resonances 
in order to avoid recurrent scattering. The calculation 
is simple for an MG dielectric. Since there are no longi- 
tudinal bulk normal modes in the regime of frequencies 
considered, coherent modes are given by the transverse 
MG polarization propagator, 



gf G {q) = £ LL G e / f (q), 



(61) 



In Eq.@U, £ LL 

factor which derives from Eq. (|60l) and 



is a Lorentz-Lorenz local field 



(62) 



is the transverse Dyson propagator of an effective 
medium with dielectric constant e. Straightforward ap- 
plication of Eqs. (|ll59l61[) yields, 



a/; 



■ Coh ,^_ 5R{e} + 2 



MG 



(63) 



For a non-absorptive medium the index of refraction sat- 
isfies n 2 = e and thus Af M °G"(r) cx (n 3 + 2n)/3. Note 
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that this dependence on n differs from that of the light 
spectrum in an effective medium (readily derivable from 
Eq. (|62]) ). J\f h 9 ht n due to the presence of a local 
field factor. The precise form of this relation is how- 
ever model-dependent. As an example, emission from an 
Onsager-Bottchcr (OB) cavity would yield Mqb 1 '^) °c 
9 ^" +1 instead [3l|, where the OB local field factor reads 
£ob = 27TT- 

From Eqs. (|6H6"3"|) it seems that long- wavelength radia- 
tive modes are given only by ^s{G e /^ (q)} in the trans- 
verse coherent emission. However, Q{ X|1 ffi ) contains also 
transverse parts of the in- free-space propagator Q. It 
can be shown that, for the long- wavelength modes of the 
effective MG dielectric, the trace of the scattering lon- 
gitudinal self-polarization propagator can be written as 



m MG 



d 3 <? [ X||(a)/pa 



(27r)3Lfc2[l + X|| ( 9 ) 
(2nY 



G 



(o) 



MG 



C LL G e l f {q)[C LL -l] 



(64) 



As a result, the total LDOS for far field emis- 
sion contains a dispersive component proportional to 
£ll,ob{£ll,ob — 1)"- in the MG and the OB models 
respectively. We can conclude that generically and for 
the far field emission in an effective medium, one local 
field factor enters the coherent spectrum while two fac- 
tors enter the total emission spectrum. Only the latter 
is well-known 13511 . 



VI. THE VACUUM ENERGY 

It has been recognized for a long time that the Casimir 
effects, the van der Waals (vdW) forces and the Lamb- 
shift share a common origin. It is customary to ascribe 
the vdW and Lamb energies to short-range interactions 
while the Casimir effects are thought of as long-range 
interactions. According to this qualitative distinction 
the vdW and Lamb energies have a microscopical ori- 
gin. It roots on the EM interactions between the dipole 
constituents of the dielectric. For instance, the sum of 
many-body interactions yields a sort of virial expansion 
for the vdW energy 36] . On the contrary, the Casimir en- 
ergy would have a macroscopical origin. It roots on the 
long-range (retarded) EM interactions between macro- 
scopic dielectrics. Typically only the long- wavelengths 
of the EM field are relevant for the Casimir energy and 
their spectrum gets discrete as the field obeys macro- 
scopic boundary conditions. Some authors refer to the 
energy of discrete modes as Casimir energy while they 
term bulk energy that of the continuous spectrum 37- 

Rather than making the distinction on the basis of the 
interaction range, an alternative terminology consists of 
referring to the ground state EM energy of the system 



when stored in the material degrees of freedom (d.o.f.) 
of the dielectric as binding energy 0, [H, [4(|. On the 
contrary, when the ground state EM energy is computed 
out of the EM vacuum fluctuations it is termed vacuum 
energy [41]. Evidently, this distinction is just conceptual 
and the total binding energy must equal the total vac- 
uum energy. In either case, the EM energy manifests in 
the shifts of bound state levels or in the shifts of the res- 
onant frequencies of optical normal modes. 

Since we deal with an infinite volume complex medium, 
surface terms and discrete modes are absent and there 
is no need for using the interaction-range-based classifi- 
cation. The computation of Q vc and the spectrum of 
vacuum fluctuations allows us to compute the total vac- 
uum energy density, J- v . On the other hand, the renor- 
malization of the single particle polarizabilities allows us 
to identify the shifts in the binding energy of atomic 
states. Without ambiguity we can term the energy shift 
of atomic levels Lamb-shift [211 |23| . Correspondingly, 
the energy density stored in the internal d.o.f. of every 
atom is the Lamb energy density, F L . It follows from the 
above reasoning that T L must be part of J- v . Nonethe- 
less, because other than internal d.o.f. are present in a 
molecular dielectric, it is expected that the binding en- 
ergies of collective (i.e . cluster) d.o.f contribute also to 
T v . In references [U EH the authors have found that, 
al leading order in pa, the Lamb-shift can be derived 
from the total variation of the vacuum energy. In turn 
this would imply that, at leading order, all the vacuum 
energy gets stored in the internal d.o.f. of the dipole con- 
stituents. However, that result is not totally conclusive 
as near field interactions and local field factors are ne- 
glected there. 

The aims of this Section are, in the first place to ob- 
tain an analytical formula for J- v which facilitates fu- 
ture investigations about its physical content. Second, 
to compute the vacuum energy of an effective medium to 
estimate in first approximation the contribution due to 
local field factors. 



The vacuum energy of a homogeneous and 
isotropic molecular dielectric 



We apply the variational method of Schwinger to the 
computation of the total vacuum energy density of a 
generic complex medium made of point dipoles. Accord- 
ing to that approach the vacuum energy emerges as an 
effective potential of the EM field. This effective po- 
tential is induced as the dipole constituents are brought 
adiabatically from infinity and are assembled until com- 
pleting the final dielectric configuration. To an infinites- 
imal variation on the dielectric properties of the medium 
corresponds a variation on the spectrum of vacuum fluc- 
tuations, J\f emls - 7 and hence a variation on the effective 
potential. The integration of the latter variations is T . 
Following [43-45], the Schwinger formula for the effective 
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action reads, 



2 IJq u > 



H 



(65) 



where the lower limit of integration can be set to zero af- 
ter regularization of divergent terms and T-L is the proper- 
time Hamiltonian operator. In our case T-L is the inverse 
of the self-polarization propagator, T-L = G 1 , in accor- 
dance to the fluctuation-dissipation theorem [Tlj . The 
Casimir energy relates to S by T v = — 9?{<S}/(VT), T 
being an asymptotic time of observation and V being 
the (infinite) volume occupied by the dielectric material. 
This way, variations on J- v are due to variations of the 
self-polarization propagator. Replacing H with G , per- 
forming the variational derivative, integrating in u and 
taking the trace, one obtains, 



5T l 



(Ilo 



.(270 7 (2tt)3 
[2Gx (q)SGl 1 (q) + G\\ (q)SG^(q)} } . (66) 



The variational derivatives can be performed w.r.t. any 
functional that characterizes the initial and final vacuum 
state. In our case, the final state is determined by G v 
which is function of G, x an d &■ The initial state is that 
in which the EM field is free and the dipoles do not in- 
teract, being decoupled from the radiation field. For the 
two-level atom model of Appendix [Bl their 'initial' polar- 
izability is a' = ^/(l — o; 2 /o; 2 ). Functional integration 
of Eq. (f6"6")) yields, 



2 

-h. 



da; 
2^ 

da; 
2^ 
da,' 
2^ 



(2tt)3 
d 3 g 
(2tt) 3 

d 3 q 
(2^)3 



In 



In 



In 



[(a'/af]} (67) 
xiGix||G||]}68) 
[Gf] 2 ]}, (69) 



where the cj, g-dependence of the integrands has been 
omitted. An expression like this was firstly obtaine d by 
Bullough and Obada for a cubic molecular crystal [4fJ]. 
The authors however failed in not considering the renor- 
malization of the single-particle polarizability. This is 
crucial in order to identify the Lamb energy. For future 
purposes we will refer to the terms which involve the 
bare and Dyson's transverse propagators in Eqs. (l68l69[) 
as Schwinger-bulk energy, 
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(70) 

Beside the in-free-space Lamb-shift Buhmann et al. 
have computed the additional shift due to the interac- 
tion of every dipole with the rest of the dielectric, £^ h 
[23l W<\ . In other words, £^ c sh amounts to the energy 
stored in the atomic d.o.f. that it would cost to take out 



of the medium one of the dipoles - note that other possi- 
ble energy costs in the removal process are not accounted 
for in £^f h ■ £sc' 1 1S made of a resonant term which only 
contributes if the dipole is in an excited state, S^es '■ 
an off-resonant term which depends on the polarizability 
of each atomic state, We just adapt the computa- 

tions in [23, 24] to our results in Sections III Bl and Mil for 
two-level atoms in the VC scenario and obtain, 



£ ™ h = ^^ fc2 [2^ c +^ni^+^}, (7i) 
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where the tp-factors are implicit functions of k. Eq.(|7ip 
equals ch(k res — ho) in Eq. (|37p except for the fact that it 
is not evaluated at k — k res but at k = k res + icT instead. 
For the same polarizability model, a\k=±iu in Eq.(|72p is 
given by Eq. (|35|l evaluated at the imaginary frequency 
ck = ±icu. The expansion of the term in Eq. (|67[) al- 
lows to verify that both the in-free-space Lamb-shift and 
p£g C sh enter that equation but with opposite sign 4 . This 
leads to an interpretation of Eqs. (|67l68l69|) slightly dif- 
ferent to that in That is, F v takes account of the 
in- free-space EM modes of Eq. (f69|) and the Lamb energy 
of Eq. ((67j) and substitutes them with the optical modes 
and/or binding energies of the coupled system of Eq . ([68|) . 

The decomposition of T v suggests two questions. 
Firstly, in view of the above preliminary results, it seems 
likely that T L is different from T v . Consequently, the 
total energy cost in the removal of one of the dipoles from 
the dielectric would be different to £^f h - This must be 
verified. Second, provided that T L 7^ F v ' , it should be 
investigated whether or not there exists any observational 
signature of such a distinction. We leave these issues for 
future work. 



B. The vacuum energy of an infinite-volume 
effective medium 



In this section we abandon the microscopical descrip- 
tion of the dielectric and concentrate on the long wave 
length modes of the vacuum energy density in an infinite- 
volume effective dielectric. The original computation car- 
ried out by Schwinger [43j was motivated by the problem 
of the sonoluminescence phenomenon on big bubbles. In 
here we will not try to address in detail this problem 
since it involves secondary issues like those of surface 
corrections [HI, etc. Our purpose is to investigate how 
the presence of local field factors, ignored in previous ap- 
proaches, modifies the usual Schwinger-bulk energy. 



4 The proof of the equality _Eo. l|67)l 
where else. 



—3J 7 will be given some- 
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The Schwinger bulk scenario restricts to low frequen- 
cies and long- wavelength modes in comparison to the typ- 
ical length scale of the dielectric microstructure. Thus, 
without much loss of generality we adopt the Maxwell- 
Garnett model of Subsection IV CI This way, the MG 
relation between e and a in Eq. (|60|) can be used and the 
Lamb-shift can be considered integrated already in a. 
Following [43j we take a large dielectric body and ignore 
the surface terms. Thus, we will refer to the resultant 
energy as effective-volume energy to distinguish it from 
the Schwinger-bulk energy and from the total vacuum 
energy which would include short wavelengths. 

Instead of using the exact formulae of Eqs. (|67l68l69[) 
it is easier here to start with the functional derivatives 
of Eq. (|66|) . For an effective medium of uniform permit- 
tivity e'(oj) = 1 + Xmg( w ) the variations on the r.h.s. of 
Eq.® stand for 



(73) 



where the explicit dependence on ui has been omitted. 
Using the formulae of Subsection IV CI we are left with 
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ft r [°° du f d'q f\ , 

x (2£ LL G e / / ( g )i 7{£z i[ Gl // ((Z )]-i } 
+ C LL Gf^{L- L l[Gf]^})}, 
where Gf, = -rrs and we have omitted the constant 

e' k z 

term — 3^5 j Jq 00 ^ In which is a reminder of the 

internal atomic binding energy of the isolated dipoles. 
The above formula presents two separable contributions, 
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The term in Eq. (|74|) is the usual Schwinger-bulk term 
equivalent to that in Eq. (|70l) for an effective medium, 
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(76) 



with n(uj) = di{ y/e(w)}. Since F^Il' nas been profusely 
studied in the literature we will not attempt to evaluate it 
further -see eg. (2jj. Nonetheless, for future discussions 



let us write T^h. as 
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bulk 



-1 



1 



do;— huWg&tp) - (77) 
By equating Eq.([76]) with Eq.([77| we define M^(w) = 



The term in Eq. (|75| comes from the combination of 
the local field factors and the effective longitudinal prop- 
agator. We will denote it by AJmg- It is not associated 
to normal modes. As a matter of fact, only longitudi- 
nal photons contribute to the LFFs in the MG model. 
Therefore, ATmg must contain the contribution of in- 
ternal resonances. We can estimate its value using a 
Lorentzian dielectric constant of band- width T and reso- 
nant frequency ui res , 



e = 1 
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(78) 



With the above formula, Eq.([75|) reads approximately, 
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(79) 



where / is an effective oscillator strength with / c 1 and 
p^ 1 is the volume-per-dipole such that J j^s = P- Note 
that the isolated contribution of LFFs would be of order 
/ instead, 



T LFF - -nh 
^MG — P n 



-^{ln[£l L ]}«^ 
(2tt) 1 1 LLU 2 



2tt 



Except for numerical pref actors, the expression for 
AJ^mg appears generic for the effective-volume energy 
density of long-wavelength modes and regardless of the 
particular microscopical dielectric model. Nonetheless, 
some comments are in order regarding the particulari- 
ties of the MG model. It was mentioned in Subsection 
IV CI that the MG model neglects recurrent scattering. 
However, at resonance scatterers overlap optically and 
recurrent scattering becomes dominant. Therefore, de- 
viations from the MG results are expected (46|. Never- 
theless, although the precise form of the LFFs changes 
at resonance, the integration over uj of Eq. (|75|) is not ex- 
pected to vary w.r.t. the result of Eq.(f?9"|) but for the 
/-dependent prefactors. 



VII. DISCUSSION ON LDOS e 



Firstly, let us comment on the distinction between 
A/^-H = 3^n 3 M and J\f llght {u) = ^n(w) for 
an effective medium. Although both spectra are com- 
puted out of G e /f and thus ascribed to the light vac- 
uum, their respective dependence on the refractive index 
and their physical meaning are quite different. We read 
from Eq.([7ijl that Afgffc oc n 3 originates from the in- 
tegration of the variations of [G^] _1 , S[G e /^ (q)}^, 1 — 
^j-<5e', weighted by the light LDOS of each intermedi- 
ate dielectric configuration of permittivity e', J\f^ ght ~ 
Q{G e /f(r, r)} e '- Because the Schwinger-bulk energy can 
be written as an integral in momenta by identifying 

" = ?S> Hulk = -eg- & the 

n 3 dependence of M^^(uj) can be fairly attributed to 
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its c 3 dependence [42|. On the contrary, J\f h 9 ht cx n 
is just the light LDOS of the final configuration of per- 
mittivity e and is simply proportional to $s{G e /* (r, r)} e . 
Alternatively, the relation M light = 3n~ 2 }f£$; (u) can 
be obtained if one starts from Fermi's Golden rule and 
restrict the quadratic vacuum fluctuations to those of the 
Dyson field. The macroscopic quantization of each field 
operator yields a factor while the remainder is the 
sum over momenta proportional to n 3 [4 71 ]. Note how- 
ever that we have defined originally the LDOSes through 
relations of the form of Eq.([T]). That is, by definition 
an LDOS is proportional to the imaginary part of some 
physically meaningful Green's function. It is in this sense 
that Nbuik' derives (i.e., is not defined) from the original 
^flight anc j likewise f or t ne LDOS of T v which derives 
from N emls - . 

Second, let us emphasize that J\f h 9 ht must not be 
identified with the total radiative LDOS. As mentioned 
in Section IIII1 both the stimulated and the sponta- 
neous emission of a point dipole are proportional to 
the quadratic fluctuations of the self-polarization field in 
\VL) emis - [TSUi]. Thus, we can write, 

Foe emls -{n\E^"{r) ■ E^"{r)\n) emls -. (80) 

The above expression together with Eqs. (|56l57[) and 
Eqs. (|22l23j) show clearly that it is not Jsf ll9htb the rel- 
evant part of the LDOS entering Femi's Golden rule in 
the descriptions of the rate of spontaneous radiative decay 
as erroneously claimed in |19| . This has been app reciated 
by a number of authors -see eg. [2(3, [H, H3, HJ HH who 
have attributed the discrepancy to the presence of LFFs 
which relate the microscopical field to the macroscopical- 
Maxwell field. We have proved microscopically that this 
is indeed the case, but rather than relating microscopical 
and macroscopical fields, LFFs in Fermi's Golden rule 
relate the self-polarization propagator to the Dyson-bulk 
propagator, they both being computed at the same mi- 
croscopical level. LFFs emerge in Fermi's Golden rule 
as a result of taking proper account of correlations. We 
found in Section IIVI that the correct statement is that 
it is only the power transferred radiatively and directly 
from the emitter into the medium that has the spectrum 
of external light. There are other contributions to the to- 
tal power which, although not transferred directly from 
the emitter into the medium, arc radiative anyways. In 
particular, the indirect coherent emission is of course ra- 
diative, although that emission is carried out by induced 
dipoles and not by the emitter itself. In the far field, 
at a distance much greater than the typical correlation 
length between the emitter and the surrounding dipoles, 
it would not be possible to distinguish between direct and 
indirect coherent radiation. 



5 flight is d enote( j by N rad in 0. Also in this context, radiative 
emission is synonym of far field emission. 



The restrictive identification in [19( of the radiative 
LDOS with J\f h 9 ht seems to have led to confusion on 
the authors of [48[ who are interested in the relation be- 
tween the refractive index and the radiative decay rate 
of an emitter embedded in a random medium. They ar- 
gue that, since in the long- wavelength limit of the effec- 
tive medium theory, q£ — > 0, the dielectric constant be- 
comes a local and isotropic function, correlation effects 
can be disregarded and the effective transverse propaga- 
tor of Eq. (jr?2"]) determines the radiative LDOS which en- 
ters Fermi's Golden rule. On the contrary, because under 
no physically acceptable condition C\ q ^^Q = 1 holds -eg. 
Cll, that argument is generically erroneous. 

A similar reasoning applies to [49| where the authors 
study the radiative decay of a dipole emitter embedded in 
a piecewise continuous dielectric. They take as LDOS of 
reference that of a homogeneous (i.e. one-piece) medium, 
~ n. However that is not the appropriate LDOS for emis- 
sion but for light propagation. The reason being that it is 
not possible to embed an emitter in a continuous medium 
and take a continuous limit such that the emitter forms 
part of the medium and the effects of LFFs disappear - 
eg. in the limit Rq — > of an Onsager cavity of radius 
R, CRq^o — £-ob 7^ 1 [IS HH- Otherwise, unphysi- 
cal divergences appear for the case that the permittiv- 
ity be a complex number. In conclusion, while J\f h9ht is 
well defined in a strictly homogeneous medium, J\f emls - is 
not. The computation of G± in [5(| for an heterogeneous 
medium is correct provided that the dielectric function 
takes proper account of the emitter embedding. How- 
ever, in such a case there are no divergences in 
and there is no need to restrict the radiative LDOS to 
as claimed by the authors. In particular, the 
statement that 3{Tr{CJ}} = 23{^} in a medium of real 
dielectric constant is in general incorrect. If properly 
computed, 2Q{G±} in an effective medium yields only 
coherent modes. 

Note however that if one were interested only in de- 
tecting the band gaps of radiative emission in an effec- 
tive medium, the detection of the gaps of J\f bght , which 
correspond to bands of total reflectivity, might suffice. 
The reason is that, for an effective medium like that of 
an MG dielectric, the total radiative spectrum is g iven 
by C 2 q (u,)A^H « ^=o(w)n(w) [!, M M HI ■ 
Therefore, unless the LFFs diverge at those frequencies 
at which the index of refraction vanishes, the gaps on 
Af h9ht (uj) coincide with those of the spectrum of total 
radiation. It is to this respect that the work of [48[ is still 
useful. In physical terms, this implies that there does not 
exist coupling of either the emitter itself or the surround- 
ing induced dipoles to propagating modes. This however 
does not preclude the non-radiative energy transfer be- 
tween the emitter and the induced dipoles. 

Finally, a comment is in order on the application of 
Eqs. (|22ll29]) to the spontaneous decay rate of a dipole 
in a statistically homogeneous complex medium. When 
the emission is stimulated, under stationary conditions 
the rate at which the dipole absorbs the energy sup- 
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plied by the external source is (ideally) balanced by the 
rate at which the dipole radiates that energy. Therefore, 
its polarizability can be considered at any time equiv- 
alent to that of the rest of dipoles in the medium and 
Eqs. (|22ll29p are strictly applicable so that N emis - oc 

Tr|3{^ yc, (r, r; cj)}|. In contrast, when the dipole de- 
cays spontaneously from an excited state, its polariz- 
ability is in general different to that of the rest of sur- 
rounding dipoles. This brakes manifestly the statistical 
translation-invariance of the medium at the emitter loca- 
tion since it behaves as an impurity. As a consequence, 
Eqs. (|22ll29[) may be applicable only in an approximate 
manner. This is at the root of the distinction between 
the virtual cavity and the real cavity scenarios [2(J [5(| ■ 



VIII. THE PROPAGATOR OF THE 
SELF-POLARIZATION FIELD IN NUMERICAL 
SIMULATIONS 

The computation of a and the measurement of all the 
components of W u can be simulated numerically using 
the so-called Coupled-Dipole-Method (CDM) [5l[. In 
that method, the transference matrix t of Subsection lll Al 
is computed for each configuration of point dipoles be- 
fore performing the average over the configuration ensem- 
ble. Let us take one of those configurations composed by 
N + l host scatterers. We decide to stimulate the dipole 
at i?o with an incident (inc) monochromatic field of fre- 
quency u) — kc, Eq. That is, Ei nc {f) = Eo if r = Ro and 
Ei n c(f) = otherwise. Generally the dipoles are not all 
equivalent and have different renormalized polarizabili- 
ties a 1 , i = 0,..,N. Therefore, the dipole moment of the 
emitter reads in terms of the incident field, 



p(R Q ) = e &°E . 



On the other hand, the total field at the emitter location 
is the sum of the incident field plus the self-polarization 
field which the dipole moment p(Ro) itself creates at Rq. 
That is, 



E(R ) = E q(Ro,Ro)-p{Ro) 

= E ~k 2 &°g(R Q ,R Q )-E Q , (82) 

where q(Rq, Rq) is the self-polarization propagator of the 
dipole for the scatterer configuration considered. The 
w-dependence is omitted hereafter. Our interest is in the 
second term on the r.h.s. of Eq. (|52"|) . It can be written 
also as the sum of all the fields which propagate from 
each of the induced dipoles to Rq , 



j 2 Jv 

k 2 &°g(R , R )-E = —Y] & 0) (Ro ~ Mi) ■ v{R\). 

i=0 



Except for the emitter, the rest of dipoles are only in- 
duced by their mutual interactions. The dipole moments 
induced can be computed formally using the t-matrix of 
the system. t(Ri,Rj) yields the dipole induced at some 
point Ri as a result of its interaction with the dipole ex- 
cited by the incident field at some point Rj. The dipole 
moment of a dipole located at Ri reads, 

, 2 _i N 

p(Ri) = -(— ) ^t(^,^)-^„ c (^-). (84) 



3=0 



Because in our case the only dipole excited externally is 
that at Ro, we get 



p(i?0 = -(^) ^(iMoMo- 



(85) 



Inserting Eq. (p5"l) into Eg. (155)) we end up with, 



1 N 

q(Ro,Ro) = ^ 1 ^Y J G ( - 0) (Ro-Ri)-t(Ri,Ro). (86) 

i=0 

Should all the dipoles be equivalent with a? = a 
Vi and density p, it is immediate to obtain 

g(Ro,Ro)) = G VC (R ,R ) as given in Eq.QZ). 



A comment is in order concerning the use of the renor- 
malized single particle polarizabilities in the above for- 
mulae. In simulations, it is recognized that the applica- 
tion of the CDM equations needs of some a priori pre- 
scription for the value of the renormalized single-particle 
polarizability in order to obtain physically acceptable val- 
ues for the scattering cross-section [52]. As a first ap- 
proximation, one can use polarizabilities which are only 
renormalized by in-free-space radiative corrections -see 
eg. [26j . A priori additional renormalizations were car- 
ried out by Draine and Goodman in 53j and Chaumet et 

(81) 

al. in f52|. Nonetheless, it must be clear by now that be- 
cause the ^-factors depend on the renormalized value of 
a and a receives radiative corrections from the (^-factors, 
the problem is one of self-consistency. In physical terms, 
this reflects the double and complementary role that ac- 
tual dipoles play in a complex medium. On the one hand, 
they polarize the vacuum. On the other hand, they renor- 
malize their own polarizabilities. 



IX. CONCLUSIONS 

We have computed the propagator of the polarization 
field in a statistically homogeneous and isotropic medium 
made of indistinguishable isotropic point dipoles. In 
Fourier-space, a relation of proportionality with Dyson's 
propagator has been found in terms of local field factors 
-Eqs. (|22l23l58l) . The stochastic kernel of its Lippmann- 
Schwinger equation is that of Eq. ([29|) . The self- 
polarization propagator determines the LDOS of the to- 
tal dipole emission, j\f emls -_ i n contrast, Dyson's propa- 
(83) gator determines only the LDOS of the power transferred 
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radiatively and directly from a dipole emitter into the 
medium. The latter is equivalent to the LDOS of the co- 
herent light which propagates through the medium from 
an external source, J\f Ught . Correspondingly, an emission 
vacuum and a light vacuum are postulated. The symme- 
try group of the latter is included in that of the former. 

An expression for the optical theorem in complex me- 
dia has been found for classical dipoles -Eq . . Formu- 
lae for the renormalized values of the resonant frequency, 
the decay rate and the bare polarizability have been de- 
rived in Eqs. (|36H55|) for a two- level atom. They agree 
with previous QM approaches. 

In the stricto sensu virtual cavity scenario, stimulated 
emission has been classified in Eqs. (|42ll45|) attending to 
its coherent /extinguished nature. Further, the coherent 
power has been decomposed in Eq. (|50p into a direct and 
an indirect component. The direct component carries 
the energy transferred radiatively and directly from the 
emitter into the medium, its spectrum is J\f h s ht and its 
field is Dyson's. The indirect component is that of the 
field radiated by the induced dipoles which is in-phase 
with the Dyson field. Only one local field factor enters 
the coherent spectrum. Correspondingly, the coherent 
dipole field is proportional to the Dyson field, being the 
constant of proportionality a local field factor, Eq. (l59"l) . 

Using Schwinger's variational method, an expression 
for the vacuum energy density has been obtained as a 
sum of different physical contributions, Eqs. (|67l68l69p . 
The vacuum energy density of an infinite- volume effective 
medium has been computed. When local field factors are 
included, additional terms proportional to the resonant 
frequency and the line- width arise -Eq. (|79p - in addition 
to the usual Schwinger-bulk term. 

The dependence on the refractive index of the density 
of different radiative modes in an effective medium has 
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been analyzed. We have found that J\f ll9ht 
n 3 , Nmq- cx (n 3 + 2n)/3, A^™ s - cx (n 5 + 4n 3 + 4n)/9, 
where the last two expressions apply to off-resonant 
Maxwell-Garnett dielectrics only since J\f Coh - and J\f emls - 
are in general model-dependent. 

It is left for future work to verify whether the conjec- 
ture of the equivalence between the Lamb-shift and the 
total variation of the vacuum energy holds. Should the 
result be negative, a physical interpretation for the dis- 
crepancy and possible observational signatures must be 
investigated. 

Also, it would be interesting to extend the present for- 
malism to more general dielectric configurations. Beyond 
dipole-aggregate models it would be necessary to develop 
a multipole expansion similar to that in (Hjj . Mul- 
tipole sources and propagators should be defined and a 
similar diagrammatic analysis to the one performed here 
carried out. 
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Appendix A: Stimulated emission from a classical 
dipole 

We adopt a classical model in which the dipole emit- 
ter is a spherical nanoparticle of permittivity e e {oS) and 
radius a, with ko«l for the frequencies of interest. 

Let us consider that the emitter is embedded in a 
generic host medium and sits at position r. It is stim- 
ulated by a stationary external field which oscillates in 
time with frequency u, Eg(r). Formally, the averaged- 
power emitted reads [56J], 



W, 



d 3 re(r-a)p»£ *(r)}, 



(Al) 



where p(r) is the polarization density induced on the 
emitter, which is proportional to E$(r) in the linear- 
small-particle approximation and is affected by self- 
polarization effects. We introduce the self-polarization 
field through the insertion of appropriate Green's func- 
tions in the above expression, 

p(r) = J d 3 r"G(r-a)eoXe 

x / A 3 r'G{r : r"-uj)-[G {0) ]- l {r"-f"-uj)-E {f"), 



so that Eq. (|Al|) reads 



Wt ot = -3{y d 3 r X ^6(r-a) J dVdV'G(r,r'; W ) 

• [&V}-\7-r<'-,oj)-Eo(r*')-EZ(?)}. (A2) 

We drop the script u hereafter unless necessary. In the 
above formulae Xe — € e{u)~ 1 is the relative electrostatic 
susceptibility of the emitter -not to be confused with the 
susceptibility of the host medium- and 

G(r)«GW(r)^[-fc 2 Xe / 6{v~a)G{v)d 3 v\ ™ (A3) 

m=0 

is the propagator which takes account of the infinite num- 
ber of self-polarization cycles which give rise to radiative 
corrections. G(r, r*) propagates virtual photons from a 
point r* inside the emitter back to another point r also 
within the emitter. All the equations above become sim- 
ple in the small particle limit, a <C fc _1 , for the electric 
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field is nearly uniform within the emitter and so are the 
polarization density and the propagator G(r, f ). The n- 
point irreducible diagrams which enter the computation 
of G can be approximated by the series of Fig0{6) in 
which the two-point correlation functions 0(r — a) ap- 
pear consecutively as factors of a product. That way the 
corresponding integrals appear untangled and the series 
becomes geometrical. Hence, Eq. (|A3[) . The underlying 
approximation reads, 

4-7T — 

/ d 3 re{r-a)G{r) ~ — a 3 G(0) 



4tt 



2<pf + v>Y> + 2 V f + c^f J -I. (A4) f) - q(0) {f _^ £ (. 



,(°) 
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the above definitions we can rewrite Eq. (|A2[) in terms of 
electrostatically renormalized operators, 



W, 



where 



I d 3 r Xe 6(r -a) J dVdV'G(r, f") 
[G't 1 ^ - r') ■ E Q {r') • ^(r)}, (A5) 



It was mentioned in Subsection III Bl that 2ipf^ and 



</?j| ' are divergent. The divergence of ip^ is cured 
by the presence of the finite radius a. Since the limit 
Lim{ / d 3 r 0(r — a)G^ t (r; a;)} = ^p-I as ka — > is con- 
ditionally convergent, it is the Heaviside function of the 



integrand that yields the finite value ^p-l 



157j . Any other 



geometry would give a different numerical value. By 
equating that result with ^a 3 ^ ' |l (leaving Sft{2<£>^} 

still free) we obtain pjj 0) = (^a 3 jfc 2 ) _1 . The net effect 
of this regularization procedure is to dress up the single 
particle susceptibility in all the in-free-space electrostatic 
corrections. This procedure is depicted in Fig H{c). That 
way we can define xjf = TT^WlXe an< ^ obtain the bare 



7JIJJ+2 

electrostatic polarizability ao(w) 



e(w) + 2' 



With 



[1 + i/sTrr ■;;}]"' 



= -t (£ -l)0(r-a) 



•« n +.••-£• + ■••yyv.=A[ l + i/3Tr{0}f 



_ WO) 

II 

1 II 

' — »=G 



(c) 



^[VV^Vf^Trli'"";} 

(d) 



Inserting the above equation into Eq. (|A5[) we obtain the 
expressions of Eqs. (|31l32[) and the renormalized polariz- 
ability of Eq. (j33"l) . 



Appendix B: Renormalization of a two-level atom 
isotropic polarizability 



The bare polarizability of a two-level atom reads, 



2u! j2 eg) p, 
Hcq Wq — to 2 



(Bl) 



where jl ® fl is the tensor product of the dipole-dipole 
transition matrix elements. The EM interactions which 
give rise to the atomic bound state have been integrated 
out and parametrized by the resonant frequency, ujq. 
This way, only coupling to radiative modes is still miss- 
ing. In a spherically symmetric state, fx <g> fx — /x 2 /3I, 



we can write a' 



with a 

3 u£— uj 2 u 



2W1 



The self- 



polarization of the atom amounts to its coupling to radia- 
tive modes. Each self-polarization cycle carries a factor 



(a) 



FIG. 4: (a) Feynman's rules for the classical regularization 
scheme of Appendix [X] (6) Diagrammatic representation of 
Eq. (|A3[) . (c) Diagrammatic representation of the dressing up 
of Xe leading to ao. Approximation symbols denote that the 
field within the emitter is taken uniform, (d) Diagrammatic 
representation of one of the self-polarization cycles which en- 
ter the series in (6). 



fc 2 Tr{o/ • G(r, r)} = 



a uj 2 



The infinite sum of cycles is a geometrical series analo- 
gous to that of Fig[4j6) for classical dipoles. It amounts 
to 1/3 times Eq. (|3"3")l with the above definitions. 
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